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The sign problem of QCD prevents standard lattice simulations to de¬ 
termine the phase diagram of strong interactions with a finite chemical 
potential directly. Complex Langevin simulations provide an alternative 
method to sample path integrals with complex weights. We report on our 
ongoing project to determine the phase diagram of QCD in the limit of 
heavy quarks (HDQCD) using Complex Langevin simulations. 


PACS numbers: 12.38.Gc, 12.38.-t, 12.38.Mh, 12.38.Aw 


1. Introduction 

The phase diagram of strong interactions is, despite various efforts, still 
largely unknown. Several states of matter are expected to be present, which 
are relevant to many phenomena such as the quark gluon plasma, neutron 
stars and the evolution of the universe after the big bang. A possible scenario 
of the QCD phase diagram is sketched in Figure|T] A theoretical prediction 
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Figure 1. A scenario of the QCD phase diagram. 


can guide heavy-ion collision experiments, which might result in the discov¬ 
ery of different states of matter. However, the sign problem leads to a path 
integral with a complex weight and thereby prevents direct determination 
using standard lattice simulation based on importance sampling. Complex 
Langevin simulations, based on stochastic quantisation, might provide a vi¬ 
able solution to sample path integrals with complex weights [1 - 16]. In the 
following we will present an update on our project to determine the phase 
diagram of heavy dense QCD (HDQCD), an approximation of QCD in the 
limit of heavy quarks, from first principles. 

2. Complex Langevin simulation 

Here, we summarise the basics of the Complex Langevin method, more 
details can be found in mumm- In analogy to the Hybrid Monte Carlo 
method, we introduce the so-called Langevin time t, which labels the evo¬ 
lution of observables and degrees of freedom in this stochastic quantisation. 
Integrating out the fermion fields leads to a path integral with a complex 
weight 

Z = J DU |det D\ e ie e " 5YM(f/) , (1) 

if the chemical potential is real and non zero, since 

[det D(n)\* = det D {—//*). (2) 

To incorporate the complex nature of the path integral in our simulations, 
we extend the gauge group from SU(3) to SL(3,C). For small step-sizes e 
the gauge links U^ x are evolved by 


Ufj, x (t + e) = R(t) U[j, x (t), 


(3) 
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where the update matrix R(t ) can be written in term of the Gell-Mann 
matrices A and stochastic Gaussian white noise r] 

R(t) = exp [zA (—e DfjS + \fk ?/)] , (4) 

where the action includes the logarithm of the determinant. Here we study 
QCD in the limit of heavy quarks, for which the fermion determinant can 
be written in terms of the (conjugate) Polyakov loops Vg and VJ 1 

det det ^1 + h e^ T Vg'j ~ det ^1 + h e _/i / T , (5) 

X 

with h = (2 k) Nt . For the gluonic part of the action we use the full Wilson 
gauge action. To avoid runaway trajectories into the non-compact exten¬ 
sion of SU(3), we apply adaptive step-size scaling [5] and adaptive gauge 
cooling mm- Too many large excursions into the imaginary directions have 
been identified to cause the Complex Langevin method to fail by converg¬ 
ing to incorrect results. It can be shown, that if the action is holomorphic 
and suitably confined in the complex extension of SU(3), Complex Langevin 
simulations are expected to converge to the correct results [9j i 10]. The log¬ 
arithm of the determinant causes poles in the derivative of the action, and 
thereby prevents the aforementioned proof to be applied. Nevertheless, re¬ 
cent work mmm has shown that especially for large quark masses this 
ambiguity will not affect Complex Langevin dynamics. We still monitor the 
distributions of the observables and the so-called unitnorm, 

unitnorm = Tr (iJU^ — , (6) 

to avoid runaway trajectories in our simulations. 

3. Numerical setup and results 

We study the phase diagram of heavy dense QCD for fixed lattice spacing 
and the simulation parameters are given in TableQ] For HDQCD the ex- 


(3 = 5.8 

H 

to 

V = 8* x N r 

k = 0.04 

n = 0.0 - 3.2 

II 

to 

1 

CO 

to 


Table 1. Summary of simulation parameters. 


pected critical chemical potential fi c (in lattice units) is related to the bare 
quark mass by 


0-c ~ m q = 


ln(2/e) = 2.53. 


( 7 ) 
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We have improved our previous results |l4j by considering larger Langevin 
trajectories, with a maximum Langevin time of 500. The interval up to 100 
Langevin time has been discarded to to remove thermalisation effects. Using 
adaptive step-sizes e we find typical values of e ~ 10~ 4 . We have determined 
the observables every 5t = 10 2 , which corresponds to approximately 100 
sweeps in between measurements. Including auto-correlation we have at 
least 2000 independent configurations for each simulation. Figure 2 shows 



0 0.0 

Figure 2. The Polyakov loop as function of T and /i. 


the expectation value of the Polyakov loop as a function of temperature 
T and the chemical potential /r. The temperature has been converted to 
physical units using the lattice spacing of a ~ 0.15 frn, which has been 
determined using the Wilson flow mm- Each black point in Figure 2 is 
the result of a dedicated simulation. The Polyakov loop shows a clear signal 
for the deconfinement transition and a transition to higher densities. At high 
densities, m q > 1, the Polyakov loop drops again. This behaviour is an 
expected lattice artefact, at which every lattice site has been filled with the 
maximum number of fermions allowed by the Pauli principle. The coloured 
surface is a cubic interpolation to the individual simulations. The resolution 
in temperature is quite limited in the fixed lattice spacing approach, since the 
temporal extent is by construction an integer. Figure 3 shows the equivalent 
plot for the susceptibility of the Polyakov loop, which directly maps out the 
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Figure 3. The Polyakov loop susceptibility as function of T and p. 


boundary of the phase diagram of HDQCD. The deconfinement transition 
appears to be quite broad, which is caused by our limited resolution for large 
temperatures and the subsequent interpolation. We find a clearer signal for 
the transition to higher densities, which almost disappears on the plotted 
scale, using the susceptibility of the fermion density 

i a In Z 

on —- 

A r T Nf dp 

4. Conclusions and Outlook 

Complex Langevin simulations provide a viable method to determine 
the phase diagram of heavy dense QCD from first principles. Further work 
includes the identification of the order of the transitions by varying the sim¬ 
ulation spatial volume and studying the Binder cumulant. Simulation at 
different lattice spacing will improve the resolution in the temporal direc¬ 
tion and allow to asses the size of lattice artefacts. The ultimate goal is 
to repeat these simulations for fully dynamical QCD [15 ; 16] and study the 
phase diagram of QCD itself. In perspective of this goal, the work here can 
be considered as blueprint for further studies and as proof of principle. 
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